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Summary 

This publication gives the complete theoretical and 
experimental results of the first phase of a joint 
government/industry /university research study on the 
vibration characteristics of twisted cantilever plates. The 
study was conducted to generate an experimental data 
base and to compare many different theoretical methods 
with each other and with the experimental results. Plates 
with aspect ratios, thickness ratios, and twist angles 
representative of current gas turbine engine blading were 
investigated. The theoretical results were generated by 
numerous finite element, shell, and beam analysis 
methods. The experimental results were obtained by 
precision machining a set of twisted plates and testing 
them at two laboratories. The second and final phase of 
the study will concern the effects of rotation. 


Introduction 

The vibration analysis of gas turbine engine blading 
has evolved from simple beam analysis to sophisticated 
finite element techniques. Although the finite element 
analysis of this blading has become commonplace, 
significant differences in the published results for various 
methods of analysis have raised some doubt concerning 
the adequacy of these methods to accurately predict the 
vibratory characteristics of highly twisted blades. 

An example of the types of differences existing in the 
literature is the predicted first bending frequency for 
pretwisted plates. The predicted dependency of frequency 
on twist angle for twisted cantilever plates spans the 
spectrum from increasing significantly to decreasing 
significantly. Likewise, major differences in predicted 
frequencies can be found for the other modes of a twisted 
cantilever plate. 

To aid in resolving such discrepancies, to validate 
analysis methods for twisted blades, and to generate an 
extensive data base on the vibrations of twisted blades, 
a joint government/industry/university research effort 
was conducted. A series of specimens were precision 
machined and carefully tested at two laboratories. 
Concurrently, analysts and researchers with known 
interests in this problem were invited to calculate the free 
vibration frequencies and mode shapes of these 


specimens by using their best available methods. Sixteen 
groups responded to the request and submitted 19 sets of 
results. 

Because of the extent of this project and the volume of 
data collected, the documentation of the research for 
phase I has been divided into four publications. 
Reference 1 compares the theoretical methods with each 
other and with the experimental results. Reference 2 
describes the experimental portion of this project and 
includes a discussion of the support compliance effect. 
Reference 3 relates the history of the problem to the 
results of this research project. Since it was felt that other 
researchers would want access to a complete set of 
experimental and theoretical results (which would be 
inappropriate in a journal article because of the 
quantity), the results are provided herein. 

The authors wish to acknowledge the following 
organizations and individuals who made substantial 
contributions to the joint research effort summarized in 
this work: University of Dayton Research Institute 
(F. Bogner); The Technical University of Denmark 
(K. Ravn-Jensen); Indian Institute of Technology- 
Madras (V. Ramamurti); Imperial College of Science and 
Technology (D.J. Ewins and B. Hillary); Institute of 
Sound and Vibration Research; University of 
Southampton (D.P. Gao and M. Petyt); General Electric 
Company/AEG (B.G. McClintick, P. Niskode, and 
M.J. Stallone); General Motors Corporation, Allison 
Gas Turbine Division (W.N. Dalton and R. Jay); United 
Technologies Corporation, Pratt & Whitney Aircraft 
Group (R.D. Aunger and J. Hurchalla); Garrett Turbine 
Engine Company (R.J. Stockton); McGraw Edison 
Company, Turbodyne Division (M. Singh); Ohio State 
University (M. Ewing, J.K. Lee, and A.J. Wang); and 
U.3. Air Force Wright Aeronautical Laboratories (B. 
Tavner). 


Symbols 

A-U method identification 

a plate length 

b plate width 

D plate flexural rigidity, D = Eh*/ 12(1 - v 2 ) 

E modulus of elasticity 


i 


h plate thickness 

/ number of finite element nodes spanwisc 
M number of spanwise nodal lines 
m number of finite element nodes chordwise 
N number of chordwise nodal lines 
n number of finite element nodes through thickness 
A frequency parameter, \ -coa 2 Sph/D 
^ twist angle (root to tip) 

v Poisson’s ratio 

p mass density 

w frequency 

Discussion 

The geometry of a twisted plate is shown in figure 1. 
Twenty geometric configurations were considered. These 
consisted of all combinations of two aspect ratios 
(i a/b = 1 and 3), two thickness ratios (b/h = 5 and 20), and 
five twist angles (^ = 0°, 15°, 30°, 45°, and 60°). 
Reference 2 also includes experimental results for 
a/b = 2. These configurations were chosen to represent 
the range of values existing in gas turbine engine blading. 
To limit the study to a manageable size and to focus on 
the important effect of pretwist, only rectangular cross 
sections were considered (hence, camber and other airfoil 
cross-section effects were ignored), and taper and 
rotational effects were not considered. The effects of 
rotation are being considered and will be discussed in a 
future publication. 

Theoretical Methods 

The 19 sets of theoretical results were derived from 
beam theory (2), shell theory (2), and finite element 
methods (15). The participants, listed in the 
acknowledgments, were asked to use their “best” 
analysis methods to obtain the first 10 frequencies and 
corresponding mode shapes for the 20 configurations 
considered, assuming complete fixity at the clamped end. 
To determine the degree of convergence, the participants 
were also asked to conduct a convergence study for one 
of the blades ( a/b~ 1, ^ = 30°, b/h = 20). Descriptions of 
their methods and mathematical models were requested. 

A basic description of each theoretical method is given 
in table I. These methods are identified by the letters A to 
S (the experimental results are identified by T and U). 
Included are the type of element, model data, and 
method of eigensolution. The models are defined by the 
notation Ixmxn , where /, m, and n represent the 
number of nodes in the spanwise, chordwise, and 
thickness directions, respectively. For plate and shell 
elements, n is always 1 and is omitted. The highlights of 
the participants’ descriptions are given here with the 


methods grouped into five categories. More detailed 
information is available from the authors of this paper. 

Finite element methods , — Finite element methods 
based on plate, shell, and solid elements were used. 

Plate elements: Methods A to G used finite elements 
based on plate theory. Methods A and B used the w ? eli- 
known finite element program NASTRAN (COSMIC 
version, April 1982 release, CTRIA2 elements). The only 
difference was that the mass matrix formulations were 
lumped in method A and coupled in method B. Method C 
used the MARC program (revision J.2-1) and the 
triangular plate element with membrane and bending 
stiffness, known as element 49. Methods D to F used 
triangular elements based on thin-plate theory. Method G 
used the quadrilateral elements of the program known as 
SAP. 

Shell elements: Methods H to K used finite elements 
based on shell theory. Method H used NASTRAN (MSC 
version, level 61, CQUAD4 elements) and included shear 
deformation effects. Method I used the quadrilateral 
thin-shell element of the MARC program. Method J used 
an eight-node quadrilateral shell element. Method K 
included thick-shell effects. 

Solid elements: Methods L to O used three- 
dimensional elements with 12, 18, 8, and 16 nodes per 
element, respectively. Except for the convergence study 
of method N, only one element w ? as used through the 
thickness in all cases. 

Other methods . — Both shell and beam theory were also 
used . 

Shed theory: Methods P and Q used direct shell theory 
analysis. Method P used the Rayleigh-Ritz method to 
solve the shallow shell theory equations. Method Q used 
the finite difference energy method to solve the 
nonshallow shell theory equations. 

Beam theory: Methods R and S used advanced beam 
theories including warping corrections. However, only 
method R included shear deformation and rotary inertia 
effects. 

Two slightly different types of “constant thickness” 
definitions were used by the participants. Those using the 
plate and shell methods defined the thickness in a 
direction normal to the local middle surface. For this 
definition, the intersection of a plane at a constant 
spanwise coordinate with the twisted blade resulted in a 
slightly nonrectangular shape (thicker at the ends). Those 
using three-dimensional models with more than one node 
through the thickness constructed their meshes such that 
this intersection resulted in a rectangle. The specimens 
were manufactured with this latter definition. To 
estimate the differences introduced into the comparisons 
as a result of this difference in thickness definition, a 
numerical study was performed on the short, thick 
configuration by using method A. Both thickness 
definitions w 7 ere used. The results showed that the 
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predicted frequencies were only significantly different for 
the a/b = 1 and <p = 6 0° case. Even in this case the 
maximum difference was less than 3 percent. The 
differences were largest for the bending modes. 

Experiment 

The specimens (fig. 2) were precision machined at the 
U.S Air Force Wright Aeronautical Laboratories from 
7075-T6 cold-drawn aluminum bar stock. The width b of 
all of the specimens was 5.08 cm (2 in). The modulus of 
elasticity was carefully measured by using the method 
described in reference 4. The specimens were then tested 
at both the U.S. Air Force Wright Aeronautical 
Laboratories and the NASA Lewis Research Center by 
holographic interferometry. A special steel holding 
fixture was designed and used to minimize the 
compliance effects of the support. Excitation was 
accomplished either with a siren (low frequency) or a 
piezoelectric crystal (high frequency). To check on the 
repeatability, at least two test runs were made by 
different technicians. The resonant frequencies were 
found to be repeatable within 0.5 percent. As previously 
mentioned, the experimental methods are described in 
detail in reference 2. 


Results 

The results of the two experimental tests and the 19 
theoretical analyses were collected and compared. The 
frequencies were put in the form of the standard plate 
nondimensional frequency parameter, X = o;a 2 'Jph/D. 
The vibratory modes of twisted plates consist of beam 
bending and torsion, chordwise bending, and extension, 
not all uncoupled from each other. However, it is useful 
to identify modes according to their predominant type of 
motion by using the nomenclature coined for beams and 
flat plates. Therefore the modes are classified as follows: 
(1) flapwise bending (IB, 2B, . . .), (2) torsion (IT, 2T, 

. . .), (3) chordwise bending (1C, 2C, . . .), (4) edgewise 
bending (IE, 2E, . . .), and (5) extension (1A, 2A, . . .). 
Although the nodal lines of normal displacement 
determined from the eigenfunctions are, in general, 
curved and do not lie exactly in the spanwise and 
chordwise directions, classifying modes by the number of 
nodal lines lying approximately in these directions is 
useful. The convention used in this study is to classify 
modes as M/N, where M and N are the number of 
normal displacement nodal lines approximately parallel 
to the spanwise and chordwise directions, respectively. 
The nodal line at the root is not counted. For the flapwise 
bending modes, M=0; for the torsional modes, M=l; 
and for the chordwise modes, M>\. For the edgewise 
bending and extensional modes, this classification is not 
appropriate. The discussion to follow uses the terms thin 


(b/h = 20), thick (b/h = 5), long (a/b = 3), and short 
(a/b= 1). 

The participants were asked to conduct convergence 
studies for one configuration. The results of these studies 
are shown in table II. Typically, three or four mesh 
densities were used and the results compared. For most of 
the methods, the difference in frequency between the 
mesh density used for the vibration study and a slightly 
coarser mesh was less than 0.2 percent for the first mode 
and less than 3 percent for the higher frequency modes. 
Examples of the finite element meshes used for the 
vibration study are shown in figure 3. 

The frequency results have been put in both tabular 
(table III) and graphical (fig. 4) form. The tables 
compare the frequencies of all methods as a function of 
twist angle for each mode of each configuration. The last 
three rows give the percent difference between the 
average of the theoretical results and the experimental 
result T. A positive number indicates that the theoretical 
average was larger than the experimental value. A dashed 
line was used when either the experimental or theoretical 
value was not available. The figures display the 
dependence of frequency on twist angle for each mode of 
each configuration. The theoretical data points are 
connected with solid lines and the experimental points are 
connected with dashed lines. The symbols represent the 
methods A to U. 

The mode shapes are also important. Five of the 
participants sent in at least partial sets of mode shapes. In 
addition, the experimental mode shapes were determined 
by holographic interferometry. The available analytical 
and experimental mode shapes total 1200 (20 specimens 
by 10 modes per specimen by 6 sets of mode shapes). 
Since it is impractical and of little value to present all 
1200 mode shapes in this publication, only those for 
the limiting twist angles (i.e., 0° and 60°) are shown in 
figure 5. 

Summary of Results 

The complete set of theoretical and experimental 
results of a joint government/industry/university 
research study of the vibration characteristics of twisted 
cantilever plates are given. The study included 19 sets of 
theoretical results obtained by using state-of-the-art 
techniques and two sets of carefully conducted 
experimental results. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, May 17, 1985 


3 



References 


1. Kielb, R.E.; Leissa, A.W.; and MacBain, J.C.: Vibrations of 
Twisted Cantilever Plates — A Comparison of Theoretical Results. 
Int. J. for Numerical Methods in Engineering, vol. 21, issue 8, 
Aug. 1985, pp. 1365-1380. 

2. MacBain, J.C.; Kielh, R.E.; and Leissa, A.W.: Vibrations of 

Twisted Cantilever Plates — Experimental Investigation. ASME 
Paper 84-GT-96, 1984. 


3. Leissa, A.W.; MacBain, J.C.; and Kielb, R.E.; Vibrations of 

Twisted Cantilever Plates — Summary of Previous and Current 
Studies. J. Sound Vibr., volT 96, no. 2, Sept. 22, 1984, pp. 
159-173. 

4. Mohr, D.Ct.: Experimental Measurements of Material Damping of 

Aluminum and Graphite/Epoxy in Free-Fall with Tuneable 
Excitation. MIT Space Systems Laboratory, Report SSL- 11 -82, 
May 1982. 


4 



TABLE I. -DESCRIPTION OF METHODS 


Elcmcni data 


Type of element 


Plate 

elements 


Shell 

elements 


Three- 

dimen- 

sional 

elements 


Other 

methods 


Triangular 

plate* 


Triangular 

plate* 

Triangular 

plate* 

Triangular 

plate 

Triangular 

plate 

Triangular 

plate 


Quadrilateral 
thin platc c 


Quadrilateral 

shell 1 ’ 

Curved quad- 
rilateral 
shelK 

Quadrilateral 

shell 


Supcrparamct- 
ric thick 
shell 


Parabolic 
conoid shell 


Quadrilateral 
thick plate 


Isoparametric 

Isoparametric 


Number 

or 

nodes 


t 

4 


Number 

of 

degrees 

or 

freedom 


t 

20 


Mass 

form 


Lumped 


Coupled 


Lumped 


Coupled 


Coupled 


Coupled 


Lumped 


Coupled 

Lumped 

Coupled 


Lumped 


Lumped 


Coupled 


Cubic for normal 
displacements; 
linear for in-plane 
displacements 


Cubic in area coor- 
dinates for normal 
displacements; 
linear for in-plane 
displacements 
Four triangles; 
cubic for normal 
displacements; 
linear for in -plane 
displacements 


Transverse shear 
included 
Double -curved, 
quadrilateral, thin- 
shel! clement 
Assumed stress hy- 
brid; modified 
Hcllinger-Rcissner 
principle; 2x2 
integration in 
plane; analytical 
through thickness 
Based on Mindlin’s 
theory including 
shear deformation 


Quadratic displace- 
ment spanwise; 
linear in other two 
directions 
Parabolic shape 
functions— in- 
plane linear 
through thickness 
Corner node brick 

14-point integration 
rule 


Model data 


a/b = 1 


/ x m x n f 


9x9 


11x11 

9x9 

11 x 9! 


11x11x2 

25 x 13x2 


Number 

of 

degrees 

of 

freedom 


360 

360 

360 

450 

252 

252 


660 

1872 


Number 

of 

elements 


128 

128 

144 

72 

72 


a /b - 3 


100 

64 


100 

72 


25x9 
25 x 9 
31x9 
13 x5 
13x5 


llxl! 
25 x 9 


1) x 1 1 x 2 
25 x 13x2 


Number 

of 

degrees 

of 

freedom 


1080 


1350 

360 

360 


550 

2592 


Number 

of 

elements 


384 

384 

480 

96 

96 


100 

192 


Method of 
eigen - 
solution 


Inverse 

power 


Inverse 

power 

inverse 

power 

Subspace 

iteration 

Sturm 

sequence 

Sturm 

sequence 


Si urm 
sequence 


Modified 

Givens 

inverse 

power 


588 


660 

1872 


100 

72 


Matrix 

iteration 


Shallow shell theory; Rayleigh-Ritz; 6x6 polynomials 

Nonshallow shell theory; finite difference energy method; 12 x 12 rectangles 

Beam theory; includes shear deformation, rotational inertia, and torsional warping; Rayleigh-Ritz, 1 0-term polynomials 
Beam theory; includes torsional warping; Galerkin technique; 15 uniform beam modes 
Experimental, Air Force Wright Aeronautical Laboratories 
Experimental, NASA Lewis Research Center 


Dynamic | 
reduc- 
tion 


Yes 

No 


Subspace 

iteration 

Simulation 

iteration 


a S\STR\x COSMIC, CTR1A2 dement 

*MARC, element 49 

c SAt> 

d MSC S' ASTR AS', CQI AD4 dement. 
e MARC. dement 4. 

m, and n represent number of nodes in tpanwisc, chord wive, and thickness directions, respectively. 
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TABLE II.— CONVERGENCE STUDY 


(1) Method G: a/b= I; b/h = 20; <p = 15° (2) Method H : a/b = 1 ; b/h = 20; ? = 60 1 ° 



Meshes studied 3 


4x4 

4x5 

5x4 

6x6 

b 10x 10 

1 

3.376 

3.376 

3.405 

3.422 

3.444 

2 

10.38 

10.55 

10.28 

10.46 

10.48 

3 

18.89 

19.13 

19.25 

19.79 

20.28 

4 

23.71 

24.81 

23.78 

25.47 

26.50 

5 

29.06 

29.68 

29.73 

30.83 

31.67 

6 

43.98 

44.42 

44.36 

44.66 

44.66 

7 

46,44 

48.56 

47.30 

51.15 

53.72 

8 

50.55 

53.71 

51.60 

56,87 

59.83 

9 

54.58 

55.80 

56.00 

57.98 

61.24 

10 

59.40 

62.64 

61.53 

66.14 

69.38 


Mode 

Meshes studied 3 

5x5 

b 10x 10 

20x20 

1 

3.2425 

3.2256 

3.2209 

2 

14.593 

14.219 

14.114 

3 

22.388 

21.927 

21.768 

4 

28.612 

27.989 

27.829 

5 

42.576 

41.152 

40.681 

6 

48.554 

45,943 

45.124 

7 

59.183 

55.854 

54.912 

8 

68.408 

63.722 

: 62.470 

9 

68.582 

67.213 

66.816 

10 

73.966 

69.362 

68.145 


(3) Method J: a/b=I; b/h = 20; y> = 30” 


Mode 




4x3 



119f : 


m 

1 


3.3365 

3.3350 

2.3350 

3.3340 

3.3336 

3.3348 

2 


14.332 

14.326 

14.339 

14.331 

14.328 

14.324 

3 

18.864 

18.702 

18.649 

18.697 

18.647 

18.628 

18.648 

4 

■ 

26.882 

26.854 

26.762 

26.749 

26.743 

26.719 

5 


34.274 

34.174 

34.242 

34.170 

34.136 

34.134 

6 

46.685 

46.381 

46.245 

46.340 

46.201 

46.133 

46.186 

7 

55.665 

55.017 

54.698 

55.002 

54.645 

54.499 


8 

63.708 

60.873 

60.133 

60,845 

60.051 

59.757 

n 

9 

66.473 

66.240 

65.864 

64.008 

63.817 

63.712 

62.882 

10 

73.497 

71.613 

71.103 

70.381 

69.504 

69.187 

69.166 


(4) Method K: a/b~ 1; b/h = 20; v> = 30” 


Mode 

Meshes studied 3 


I x 1 

2x2 

b 3 x 3 

■D 

n 

3.19297 

3.07876 

3.06964 

3.06597 

mm 

8.25274 

13.4014 

13.3656 

13.3367 

3 

22.1768 

18.0637 

17.3390 

17.2032 


35.5662 

25.6697 

25.7292 

25.5695 


38.8702 

32.8606 

32.4942 

32. 1968 


62.0406 

45.2402 

44.3371 

43.9086 


100.143 

55.2076 

52.1282 

50.3748 

8 

124.819 

63.6912 

57.3779 

55.2743 

9 

174.390 

70.3630 

63.8795 

61.2159 

10 

240.252 

1 81.5775 

69.8640 

66.4574 


(5) Method L: a/b = 1 ; b/h = 20; y> = 30” 


Mode 

Meshes studied 3 

m 

m 

4x8 

6x 12 

I 

3.398 

3.381 

3.379 

3.374 

2 

13.98 

14.18 

14.26 

14.31 

3 

19.21 

18.49 

18.40 

18.37 

4 

24.33 

25.21 

25.78 

26.27 

5 

33.96 

33.94 

34.03 

34.09 

6 

46.36 

46.95 

46.94 

46.88 

7 

50.30 

52.21 

52.88 

53.51 

8 

57.80 

58.39 

' 59.57 

60.87 

9 

62.60 

63.54 

61.16 

60.67 

10 

65.10 

70.20 

68.48 

68.51 

11 

86.22 

83.74 

86.04 

88.86 

12 

76.94 

92.09 

91.98 

93.09 
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^Number of elements in spanwisc, ehordwise, and thickness directions, respectively. 
1 Mesh used in vibration study. 

























TABLE II.— Concluded 


(6) Method N: a/b = 1 ; b/h = 20; *? = 30° ( 7 ) Method O; a/b = 2; b/h = 20; - 30° 


Mode 

Meshes studied 3 


4x2 

6x3 

8x4 

10x5 

b 12x6 

1 

3.812 

3.744 

3.710 

3.710 

3.710 

2 

20.89 

20.48 

20.36 

20.29 

20.26 

3 

22.36 

21.07 

20.68 

20.53 

20.45 

4 

66.20 

59.59 

58.36 

57.78 

57.51 

5 

63.08 

60.37 

59.61 

59.29 

59.10 

6 

78.62 

71.13 

69.52 

68.76 

68.66 

7 


105.04 


101.10 

100.56 

8 


107.94 

104.98 

103.77 

103.16 

9 


136.77 


127.96 

126.58 

10 



140.45 

137.89 

136.82 


Mode 

Meshes studied 3 


5x5x1 

b 10x 10 x 1 

15x15x1 

20 x 20 X 1 

10x 10x2 

1 

3.37 

3.35 

3.35 

3.35 

3.38 

2 

14.14 

14.23 

14.25 

14.26 

14.28 

3 

18.54 

18.27 

18.24 

18.24 

18.49 

4 


25.91 

26.21 

26.31 

26.12 

5 

33.29 

33.37 

33.54 

33.58 

33.67 

6 

45.71 

45.76 

45.80 

45.80 

45.95 

7 

53.38 

52.52 

52.81 

52.94 

53.22 

8 

56.72 

58.31 

58.22 

58.22 

59.09 

9 

60.46 

58.67 

59.60 

59.96 

59.64 

10 

70.29 

65.96 

66.11 

66.19 

67.10 


(8) Method P: a/b = 1 ; b/h = 20; v = 30° (9) Method Q: a/b = 1 ; b/h = 20; = 30° 


Mode 

Meshes studied 0 


Mode 

Meshes studied d 


4x4 

b 6x6 

6x8 

5x8 

4x4 

8x8 

b 12x 12 

16 x 16 

1 

3.2812 

3.2653 

3.2645 

3.2659 


1 

3.37329 

3.39250 

3.39722 

3.39929 

2 

16.444 

16.079 

16.070 

16.077 


2 

13.37254 

14.18996 

14.32722 

14.37397 

3 

19.888 

19.409 

19.404 

19.410 


3 

18.52758 

18.64300 

18.68638 

18.70206 

4 

31.467 

28.120 

27.937 

27.942 


4 

30.18277 

27.66221 

27.38769 

27.30346 

5 

38.357 

37.347 

37.318 

37.352 


5 

35.10403 

34.89516 

34.89738 

34.89309 

6 

50.347 

49.037 

49.949 

49.098 


6 

48.10618 

46.67379 

46.51712 

46.47342 

7 






7 

62.90461 

57.08569 

56.64977 

56.51397 

8 

76.455 

67.042 

65.144 

65.207 


8 

68.50919 

61.70973 

61.16853 

60.99656 

9 






9 

81.25857 

68.34251 

65.81152 

64.96341 

10 






10 

104.05667 

72.75561 

71.48077 

71.15059 


a Number or elements in spanwise, chordwisc, and thickness directions, respectively. 
b Mesh used in vibration study. 

c Effective number of terms in polynomial expansions in spanwise and chordwisc directions, respectively. Symmetric and antisymmetric analyses were made with half the number or polynomial 
terms in chordwisc direction. 

^Number of dcscritized rectangles in the spanwise and chordwise directions. 
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TABLE III.— FREQUENCY PARAMETERS 


(1) First bending mode: 0/0; a/b = 1; b/h = 20. 



Experiment 


Difference between 
average of the- 
oretical results 
and experimental 
result T, percent 



9 3.366 3.381 3.401 


3 2.965 2.842 


\9 2.965 2.853 

» 9.92 11.55 

10.44 12.75 










































TABLE III. — Continued. 


(2) First torsion mode: 1/0; a/b — 1; b/h = 20. 






















TABLE III. —Continued. 


(3) Second bending mode: 0/1; a/b = 1; b/h = 20. 


Type of data 

Method 

Twist angle, deg 

0 

15 

30 

45 

60 

Finite element 

A 

B 

C 

D 

E 

F 

G 

H 

1 

J 

K 

L 

M 

N 

0 

21.420 

20.750 

21.610 

20.660 

20.370 

21.090 

20.970 

21.270 

21.290 

21.200 

21.200 

20.810 

20.750 

20.150 

20.920 

19.970 

19.620 

18.940 

18.490 

19.110 

18.180 

17.230 

17.190 

18.270 

18.700 

19.180 

18.650 

17.340 

18.490 

16.730 

16.390 

16.910 

15.980 

15.580 

14.640 

14.360 

14.810 

13.880 

13.630 

13.600 

12.110 

14.220 

15.130 

14.400 

11.570 

14.110 

20.280 

20.560 

15.380 

16.410 

20.380 

20.080 

20.210 

16.470 

14.100 

16.270 

20.730 

22.240 

20.080 

21.710 

18.270 

20.150 

15.970 

18.070 

13.740 

16.020 

Shell 

p 

Q 

21.310 

21.290 

20,850 

20.570 

19.410 

18.690 

17.180 

16.390 

14.460 

14.210 

Beam 

R 

S 

20.730 

21.020 

19.980 

20.430 

18.210 

18.900 

16.190 

16.980 

14.310 

15.230 

Experiment 

T 

U 

19.870 

19.950 

18.000 

18.060 

16.920 

16.930 

14.540 

14.530 

12.390 

12.450 

Difference between 
average of the- 
oretical results 
and experimental 
result T, percent 

All 

5.87 

11.80 

8.65 

10.87 

12.35 

Finite element 
methods (A to O) 

5.90 

11.73 

8.25 

10.18 = 

11.60 

Selected methods 
(D and H to N) 

5.48 

10.95 

8.05 

8.36 

10.63 




























TABLE III. — Continued. 


(5) Second torsion mode: 1/1; a/b = \\ b/h -20. 


Type of data 

Method 


Twist angle. 

P, deg 




0 

D 



60 

Finite element 

A 

31.240 

32.670 

36.260 

40.590 

44.640 


B 

30.180 

31.540 

34.920 

38.940 

42.600 


C 

31.590 

33.030 

36.680 

41.090 

45.310 


D 

30.440 

31.940 

35.670 

40.130 

44.170 


E 

F 

30.080 

30.070 

32.350 

37.070 

36.610 

42.240 

41.780 

42.680 

43.800 


G 

30.550 

31.670 

34.520 

37.900 

41.650 


H 

I 

30.500 

30.970 

31.690 

34.640 

35.150 

38.120 

41.150 

41.450 


J 

31.180 

31.330 

34.170 

37.480 

40.320 


K 

30.480 

31.110 

32.490 

33.230 

32.350 


L 

M 

29.840 

30.990 

33.940 

37.570 

40.940 


N 

29.750 

30.790 

33.360 

36.470 

39.200 


0 

31.580 

32.700 

35.520 

38.980 

42.200 

Shell 

p 

30.980 

32.570 

37.350 

46.450 

62.180 


Q 

30.980 

32.110 

34.900 

38.170 

41.000 

Beam 

R 

30.910 

38.870 

55.110 

73.770 

93.420 


S 

30.870 

38.710 

54.760 

73.960 

93.430 

Experiment 

T 

28.410 

27.990 

31.330 

34.560 

37.030 


U 

28.510 

28.090 

31.290 

34.510 

37.100 

Difference between 

All 

7.39 

14.55 

16.22 

20.27 

23.60 

dvcl d.g c Vi IIIC 







oretical results 
and experimental 
result T, percent 

Finite element 
methods (A to O) 

7.17 

12.03 

10.67 

10.95 

10.99 








Selected methods 
(D and H to N) 

6.70 

10.60 

8.40 

7.01 

7.29 





























TABLE III.— Continued. 


(7) Chordwise bending mode: 2/1 ; alb = 1 ; b/h = 20. 


Type of data 

Method 


Twist angle, 

deg 




o 

15 

30 

45 

60 

Finite element 

A 


55.720 

57.710 

58.220 

58.770 


B 

51.890 

52.770 

54.500 

55.450 

56.640 


C 

55.590 

56.550 

58.600 

59.030 

59.490 


D 

53.380 

54.330 

55.940 

63.190 

67.630 


E 

54.810 

55.850 

57.150 

57.500 

58.410 


F 

53.520 

53.310 

55.760 

57.900 

58.400 


G 

53.050 

53.720 

54.920 

55.710 

57.000 


H 

52.920 

53.810 

55.660 

55.960 

55.350 


I 

54.190 


57.060 

— 

56.350 


J 

51.960 

52.890 

54.650 

54.840 

54.840 


K 

53.640 

53.650 

52.130 

47.090 

45.350 


L 

49.140 

50.140 

52.210 

52.610 

52.070 


M 







N 

50.370 

51.880 

52.510 

52.040 

51.770 


O 

55.030 

55.880 

57.850 

58.660 

58.110 

Shell 

p 

55.000 

56.120 

59.050 

63.460 

77.370 


Q 

54.620 

55.380 

56.650 

56.230 

55.910 

Beam 

R 

S 



— 



Experiment 

T 

U 

49.780 

47.540 

50.030 

55.010 

57.030 







Difference between 

All 

6.56 

12.18 

10.30 

2.68 

1.19 

average nf the- 







ur VI U^Vr V llilr 

oretical results 

Finite element 

6.14 

11.77 

9.82 

1.79 

-1.04 

and experimental 

methods (A to 0) 






result T, percent 








Selected methods 

4.69 

9.93 

7.88 

-1.33 

-4.13 


(D and H to N) 







14 





















TABLE III. — Continued. 

(9) Third bending mode: 0/2; a/b= I; b/h = 20, 







































































Type of data 


Method 


Finite element 


Shell 



Difference between 
average of the- 
oretical results 
and experimental 
result T, percent 


A 

B 

C 

D 

E 

F 

G 

H 

I 

J 

K 

L 

M 

N 

O 

P 

Q 

R 

S 

T 

U 

Ail 


Finite element 
methods (A to O) 

Selected methods 
(D and H to N) 


. — Continued. 


: 0/0; a/b = 1; b/h = 5 


















Type of data 


Method 




0 

Finite element 

A 

7.947 


B 

7.807 


C 

8.590 


D 

8.272 


E 

8.537 


F 

8.481 


G 

8.455 


H 

7.541 


I 

8.510 


J 



K 

7.354 


L 

7.211 


M 

7.556 


N 

7.127 


O 

7.597 

Shell 

p 



Q 


Beam 

R 

7.734 


S 

8.306 

Experiment 

T 

6.763 


U 

6.783 

Difference between 
average of the- 

All 

14.81 



oretical results 

Finite element 

14.69 

and experimental 
result T, percent 

methods (A to O) 





Selected methods 
(D and H to N) 

11.63 


a/b = 1; b/h = 5. 


Twist angle, 

deg 


15 

30 

45 

60 

7.859 

8.269 

8.414 

8.863 

7.722 

8.128 

8.274 

8.722 

8.667 

8.872 

9.143 

9.411 

8.321 

8.481 

8.698 

8.911 

8.612 

8.912 

8.760 

9.286 

9.810 

9.380 



8.488 

8.609 

8.823 

9.073 

7.621 

7.829 

8.630 

8.100 

8.361 

8.906 



7.326 

7.196 

6.898 

6.398 

7.283 

7.484 

7.781 

8.124 

7.734 

8.154 

8.635 

9.093 

7.191 

7.370 

7.624 

7.920 

7.683 

7.917 

8.247 

8.613 






8.593 


■ 




7.948 

8.532 

9.404 

BHPfilj 

8.524 

9.138 

10.010 


6.712 

6.903 

7.236 

1 

6.736 

6.939 

7.257 

■ 

15.33 

16.70 

15.11 

mi 

14.77 

15.68 

13.10 


11.44 

12.37 

9.05 

10.16 
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TABLE III. — Continued. 


(13) First edgewise bending mode: a/b= 1; b/h — 5. 


Type of data 

Method 


Twist angle, 

, deg 




0 

15 

30 

45 

60 

Finite dement 

A 

11.110 

10.960 

10.560 

9.980 

9.318 


B 

11.110 

10.950 

10.540 

9.943 

9.269 


C 

11.170 

1 1 .070 

10.800 

10.380 

9.855 


D 

11.130 

10.990 

10.600 

10.020 

9.340 


E 

11.570 

11.460 

11.010 

10.330 

9.800 


p 

11.570 


11.070 


9.820 






G 

10,930 

10.800 

10.420 

9.721 

8.604 


H 

I 

J 

K 

10.930 

11.120 

10.750 

10.250 

10.670 

9.577 

8.825 

9.509 





10.950 

10.620 

9.802 

8.803 

7.894 


L 

11.340 

1 1 .220 

10.880 

10.370 

9.778 


M 

10.930 

10.670 

10.060 

9.364 

8.715 


N 

10.940 

10.790 

10.370 

9.753 

9.019 


0 

10.950 

10.860 

10.580 

10.170 

9.679 

Shell 

p 






Q 



10.410 








Beam 

R 

10.730 

10.630 

10.350 

9.939 

9.458 


S 

16.770 

16.070 

14.760 

13.410 

12.250 

Experiment 

T 

U 

4.069 

5.487 

6.044 

6.934 

8.823 







Difference between 
average of the- 

All 

64.47 

51.33 

43.89 

31.52 

6.59 

Finite element 

63.42 

49.79 

42.68 

29.73 

4.56 

oretical results 

and experimental 

methods (A to O) 






result T, percent 

Selected methods 
(D and H to N) 

63.17 

49.38 

41.75 

28.13 

2.09 




TABLE III.— Continued. 


(14) Second bending mode: 0/1; alb = 1; b/h = 5. 


Type of data 

Method 

Twist angle, deg 

0 

15 

30 

45 

60 

Finite element 

A 

B 

C 

D 

E 

19.100 

18.540 

21.600 

19.630 

20.330 

1 fV\f\ 

19.110 

18.550 

21.540 

19.590 

20.220 

20.440 

20.870 

18.530 

19.090 

18.560 

21.350 

19.470 

20.220 

20.440 

20.660 

18.630 

20.910 

18.980 

18.480 

21.040 

19.230 

20.070 

18.720 

18.250 

20.620 

18.820 

19.920 

19.990 

20.030 

18.460 

19.970 


r 

G 

H 

1 

J 

K 

L 

M 

N 

0 

/ I . vUv 

18.480 

21.290 

17.710 

17.180 

18.400 

16.500 

18.670 

20.410 

18.640 



17.380 

17.170 

18.600 

16.470 

18.670 

16.460 

17.140 

18.980 

16.360 

18.680 

15.130 

17.040 

19.280 

16.180 

18.640 

13.670 

16.880 

19.370 

15.930 

18.530 

Shell 

P 

Q 






— 


20.930 






Beam 

R 

S 

17.520 

21.020 

17.590 

21.960 

17.770 

23.960 

18.000 

26.460 

18.230 

29.080 

Experiment 

T 

U 

13.900 

13.040 

13.150 

13.880 

12.430 





Difference between 
average of the- 
oretical results 
and experimental 
result T, percent 

All 

27.34 

31.77 

32.18 

27.38 

35.11 

Finite element 
methods (A to O) 

27.26 

31.41 

31.04 

25.35 

32.85 

Selected methods 
(D and H to N) 

24.68 

27.38 

28.06 

21.06 

29.32 
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TABLE III.— Continued. 


(15) Chordwise bending mode: 2/0; a/b= 1; b/h-5. 


Type of data 

Method 

Twist angle, <p> deg 

0 

15 

30 

45 

60 

Finite element 

A 

B 

C 

D 

E 

F 

G 

H 

I 

J 

K 

L 

M 

N 

O 

27.010 

26.050 

31.590 

28.650 

30.780 

30.550 

26.430 

25.470 

31.380 

28.360 

30.250 

29.420 

30.250 

25.110 

26.540 

25.540 
30.790 

27.770 
29.500 
29.650 
29.630 

24.770 
29.880 

25.570 

24.620 

29.940 

26.910 

28.600 

29.200 

28.540 

24.230 

25.300 

24.340 

28.940 

25.870 

27.550 

27.180 

27.760 

23.530 

27.330 

25.230 

30.970 



24.060 

23.210 

25.120 

22.200 

25.100 

23.800 

23.130 

24.930 

22.080 

25.030 

23.090 

22.880 

24.480 

21.690 

24.810 

22.120 

22.530 

24.000 

21.120 

24.440 

21.220 

22.110 

23.690 

20.470 

23.940 

Shell 

p 

Q 








29.910 






Beam 

R 

S 

28.130 

29.650 

28.670 

30.190 

30.190 

31.740 

32.530 

34.220 

35.460 

37.220 

Experiment 

T 

U 

21.520 













Difference between 
average of the- 
oretical results 

AH 

20.94 









Finite element 
methods (A to 0) 

20.19 





and experimental 





i C5uu i , pcitcni 

Selected methods 
(D and H to N) 

16.05 




















TABLE III, — Continued. 


(17) Extensional mode: a/b- 1; b/h = 5. 


Type of data 

Method 

Twist angle, deg 



0 

15 

Hi 

45 

60 

Finite element 

A 

— 

26.010 


25.140 

24.440 


B 

26.120 

26.020 

25.730 

25.160 

24.480 


C 

26.190 

26.110 

25.820 

25.330 

24.590 


D 

26.220 

26.110 

25.780 

25.240 

24.510 


E 

26.250 

26.050 

25.750 

25.080 

24.100 


F 

26.250 


25.380 


23.580 






G 


26.310 

25.700 

25.070 

24.770 


H 

26.150 

26.050 

25.730 

25.210 

24.500 


I 

26.500 


26.060 


24.830 


j 




K 

26.140 

25.700 

24.570 

23.010 

20.950 


L 

27.480 

27.390 

27.130 

26.690 

26.090 


M 

26.130 

26.120 

25.950 

25.570 

25.120 


N 

26.180 

26.080 

25.780 

25.300 

24.660 


0 

26.240 

26.180 

26.000 

25.720 

25.350 

Shell 

p 






o 



25.710 









Beam 

R 






s 










' 



Experiment 

T 






U 












Difference between 

All 






average of the- 
oretical results 






Finite element 
methods (A to 0) 






and experimental 
result T, percent 







Selected methods 








(D and H to N) 









































TABLE III. — Continued. 

(19) Chordwise bending mode: 2/1 ; a/b= I; b/h- 5. 


















































































TABLE III. — Continued. 


(23) First bending mode: 0/0; a/b = 3; b/h = 20. 



3 < 


a 




























TABLE III. — Continued. 


(24) First torsion mode: 1/0; a/b = 3; b/h- 20. 


Type of data 

Method 


Twist angle, <p 

deg 


0 

15 

30 

45 

60 

Finite element 

A 

21.120 

22.470 

26.040 

30.950 

36.500 


B 

20.800 

22.130 

25.640 

30.480 

35.950 


C 

21.190 

— 

26.140 

— 

36.650 


D 

21.080 

22.450 

26.030 

30.990 

36.590 


E 

21.240 

24.450 

32.090 

41.120 

51.890 


F 

21.250 


31.520 


50.220 


G 

21.690 

23.850 

27.000 

30.850 

34.370 


H 

20.540 


— 


34.240 


I 

21.140 


25.550 


33.490 


j 

K 

20.980 

22.080 

24.980 

28.910 

33.220 


L 

20.410 



28.780 



lYl 

N 

20.850 

22.010 

25.100 

29.440 

34.430 


O 

21.200 

22.200 

25.400 

29.700 

34.700 

Shell 

p 

21.150 

22.380 

26.330 

34.290 

50.770 


Q 

21.220 

22.280 

25.150 

29.200 

33.900 

Beam 

R 

21.040 

22.260 

25.550 

30.210 

35.670 


S 

21.040 

22.240 

25.490 

30.300 

35.730 

Experiment 

T 

19.580 

20.680 

23.530 

28.270 

32.290 


U 


20.790 

23.610 

28.310 

32.280 

Difference between 

All 

7.01 

8.36 

11.32 

9.31 

15.07 

average of the- 







oretical results 

Finite element 

6.93 

8.92 

12.41 

9.53 

14.32 

and experimental 

methods (A to O) 






result T, percent 








Selected methods 

6.02 

6.76 

7.42 

4.27 

6.12 


(D and H to N) 









TABLE III. — Continued. 


(25) Second bending mode: 0/1; a/b~ 3; b/h = 20. 


Type of data 

Method 

Twist angle, deg 

0 

15 

30 

45 

60 

Finite element 

A 

B 

C 

D 

E 

F 

21.380 

21.310 

21.410 

20.580 

21.580 
22.140 
21.120 
21.590 
21.360 

20.720 

20.660 

19.100 

19.040 

19.090 

18.330 

19.390 

19.060 

16.460 

17.150 

17.110 

15.320 

15.280 

15.280 

14.560 
15.680 
15.690 
11.620 
15.220 

17.560 

19.930 

21.070 

20.690 

18.790 

16.430 

17.370 


G 

H 

I 

J 

14.110 


20.220 





K 

L 

M 

N 

0 

21.380 

21.320 

20.650 

18.950 

16.900 

17.000 

14.920 



15.140 

16.300 

21.300 

22.900 

20.640 

22.300 

19.010 

20.400 

17.040 

18.300 

Shell 

p 

Q 

21.290 

21.360 

20.840 

20.670 

19.310 

19.000 

17.090 

16.990 

14.360 

15.070 

Beam 

R 

S 

20.990 

21.020 

20.320 

20.430 

18.680 

18.900 

16.720 

16.980 

14.860 

15.230 

Experiment 

T 

U 

20.760 

19.680 

19.760 

18.030 

18.120 

f — • — 

15.990 

16.140 

14.090 

14.100 


Difference between 
average of the- 
oretical results 
and experimental 

rpQtilf T nprrpnt 

All 

3.08 

4.43 

5.09 

5.16 

6.88 

Finite element 
methods (A to O) 

3.40 

4.49 

5.13 

4.95 

7.39 

1 VOUll 1 j UVi tvil l 

Selected methods 
(D and H to N) 

2.33 

3.56 

5.74 

5.06 

8.98 




TABLE IIL— Continued. 


(26) Third bending mode: 0/2; afb = 3; b/h = 20. 


Type of data 

Method 

Twist angle, sp» deg 



0 

15 

30 

45 

60 

Finite element 

A 

60.150 

58.440 

56.210 

53.840 

51.280 


B 

59.790 

58.130 

55.900 

53.340 

50.990 


C 

60.300 


56.350 


51.380 


D 

57.860 

56.450 

54.240 

51.840 

49.240 


E 

62.350 

61.000 

58.640 

55.450 

53.250 


F 

64.370 


57.970 


53.260 


VJ 

H 

62.200 





52.510 


I 

60.020 


59.150 


58.520 


J 

K 

60.100 

58.070 

55.480 

52.620 

49.430 


L 

59.840 


— 

53.410 



M 

N 

60.240 

58.370 

56.110 

53.710 

51.080 


0 

62.460 

61.470 

59.790 

57.700 

55.130 

Shell 

p 

60.100 

58.320 

56.610 

54.380 

50.340 


Q 

60.120 

58.150 

55.950 

53.580 

50.950 

Beam 

R 

58.640 

57.040 

54.910 

52.620 

50.120 


S 

58.860 

58.020 

56.780 

55.820 

55.330 

Experiment 

T 

58.650 

51.610 

46.290 

47.310 

47.680 


U 


55.840 

40.160 

49.900 

47.540 

Difference between 

All 

3.00 

11.77 

18.39 

12.43 

8.64 

average of the- 







oretical results 

Finite element 

3.55 

12.30 

18.77 

12.37 

8.96 

and experimental 

methods (A to O) 






result T, percent 








Selected methods 

2.32 

10.45 

17.70 

10.56 

8.58 


(D and H to N) 









TABLE III. — Continued. 


(27) First edgewise bending mode: a/b = 3; b/h -20. 


Type of data 

Method 


Twist angle, 

ip, deg 




0 

15 

30 

45 

60 

Finite element 

A 

63.470 

66.720 

73.000 

80.610 

87.910 


B 

63.470 

66.660 

72.920 

80.490 

87.700 


c 

63.540 


73.100 


88.080 






D 

63.690 

66.620 

72.800 

80.270 

86.980 


E 

67.400 

73.460 

83.570 

91.660 

99.070 


p 

67.460 


82.900 


98.400 


G 




H 

62.470 




87.930 



71.680 


I 

j 

62.770 

83.140 





K 

62.440 

65.680 

71.340 

77.340 

83.710 


L 

62.900 



80.280 



M 

N 





62.710 

66.000 

72.090 

79.440 

86.800 


O 

64.820 

67.340 

72.980 

80.340 

88.150 

Shell 

p 

62.290 

66.160 

74.600 

89.410 

115.600 


Q 

61.820 

65.290 

71.400 

78.820 

86.020 

Beam 

R 

61.920 

65.100 

71.380 

79.060 

86.380 


S 

67.080 

70.130 

77.460 

87.920 

99.890 

Experiment 

T 


57.820 

57.980 

64.870 

73.470 


U 


60.480 

63.180 

70.810 

76.740 

Difference between 
average of the- 

All 


13.95 

22.04 

21.02 

19.31 

Finite dement 
methods (A to O) 


14.34 

22.32 

20.21 

17.35 

oretical results 

and experimental 
result T, percent 

Selected methods 


12.53 

19.45 

18.23 

14.28 




(D and H to N) 








TABLE III. — Continued. 


(28) Second torsion mode: 1/1; a/b-^\ b/h = 20. 


Type of data 

Method 

Twist angle, deg 



0 

15 

30 

45 

60 

Finite element 

A 

66.420 

69.690 

78.500 

91.090 

105.600 


B 

65.380 

68.590 

77.300 

89.630 

103.900 


C 

66.700 


78.890 

— 

106.100 


D 

66.050 

69.370 

78.320 

90.980 

105.600 


E 

73.460 

76.830 

88.960 

109.200 

130.900 


F 

73.450 


87.610 

— 

129.390 


o 

H 

65.220 





100.600 


I 

i 

66.370 


77.870 


99.330 


J 

K 

65.890 

68.520 

75.570 

85.300 

96.090 


L 

64.080 



85.380 



IVI 

N 

65.490 

68.350 

76.270 

87.490 

100.600 


O 

67.000 

69.240 

77.320 

88.220 

101.100 

Shell 

p 

66.420 

69.430 

79.400 

100.100 

144.100 


Q 

66.680 

69.320 

76.570 

87.050 

99.410 

Beam 

R 

66.050 

69.570 

79.130 

92.750 

108.850 


S 

66.040 

69.510 

78.950 

93.000 

109.000 

Experiment 

T 

62.430 

64.140 

71.610 

82.700 

94.020 


U 


64.300 

71.700 

82.760 

93.920 

Difference between 

All 

6.71 

8.18 

9.73 

9.80 

14.04 

average of the- 







oretical results 

Finite element 

7.00 

8.48 

10.11 

9.03 

12.30 

and experimental 

methods (A to O) 






result T, percent 








Selected methods 

4.71 

6.70 

7.01 

5.26 

6.40 


(D and H to N) 
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TABLE III.— Continued. 


(29) Third torsion mode: 1/2; alb = 3; b/h- 20. 


Type of data 

Method 


Twist angle, <p 

deg 




0 

15 

30 

45 

60 

Finite element 

A 

120.400 

124.000 

134.100 

149.100 

167.200 


B 

118.300 

121.900 

131.800 

146.500 

164.100 


C 

121.000 


134.900 


168.100 


D 

118.800 

122.600 

132,900 

148.200 

166.500 


E 

122.000 

130.700 

140.200 

156.400 

173.900 


F 

G 

124.000 


137.500 


171.200 


H 

120.500 

— 



162.700 


I 

j 

120.000 


135.100 

— 

164.500 


K 

119.300 

122.400 

130.600 

142.400 

155.700 


L 

M 

115.800 



141.400 



N 

118.300 

121.500 

130.500 

143.900 

159.800 


0 

129.000 

128.600 

134.400 

147.200 

162.800 

Shell 

p 

120.100 

123.600 

135.700 

161.900 

222.000 


Q 

120.800 

123.800 

132.400 

145.200 

160.700 

Beam 

R 

119.200 

124.700 

139.700 

161.400 

187.300 


S 

119.200 

124.600 

139.400 

161.700 

187.500 

Experiment 

T 

113.500 

114.700 

123.400 

136.600 

151.200 


U 


115.000 

123.800 

136.900 

151.000 

Difference between 

avpriiop nf thp- 

All 

5.75 

7.80 

8.55 

9.20 

11.89 

dvciagc Ul llic 







oretical results 

Finite element 

5.90 

7.89 

8.05 

7.00 

8.44 

and experimental 

methods (A to 0) 






result T, percent 








Selected methods 

4.45 

6.11 

6.71 

5.12 

6.57 


(D and H to N) 









TABLE III. — Continued. 


(30) Fourth bending mode: 0/3; afb = 3; b/h = 20. 


Type of data 

Method 


Twist angle, 

deg 


0 

15 

30 

45 

60 

Finite element 

A 

118.500 

118.100 

117.300 

116.300 

116.200 


B 

117.400 

117.100 

116.200 

115.300 

115.200 


C 

119.000 


117.800 

— 

116.700 


D 

113.900 

113.600 

112.800 

112.100 

112.600 


E 

128.700 

124.000 

124.000 

122.700 



F 

r; 

128.700 


120.630 


122.650 


H 

127.100 





123.100 


I 

T 

118.000 


123.300 

— 

129.500 


K 

119.500 

118.900 

117.000 

114.500 

112.100 


L 

117.700 



115.700 



N 

120.400 

120.000 

119.000 

117.500 

116.500 


O 

122.800 

125.800 

127.600 

126.300 

125.200 

Shell 

p 

125.600 

126.200 

128.600 

133.900 

146.300 


Q 

118.600 

118.300 

117.400 

116.200 

115.800 

Beam 

R 

114.600 

114.300 

113.600 

112.600 

113.300 


c 







o 






Experiment 

T 

115.200 

112.700 

112.300 

109.100 

108.100 


U 


113.300 

112.800 

110.200 

108.600 

Difference between 

All 

4.56 

5.79 

6.13 

7.90 

10,21 

average of the- 







oretical results 

Finite element 

4.77 

5.80 

6.07 

7.19 

9.14 

and experimental 

methods (A to O) 






result T, percent 








Selected methods 

3.54 

4.09 

4.85 

5.09 

8.98 


(D and H to N) 






































TABLE III.— Continued. 


(32) Chordwise bending mode: 2/0; afb- 3; b/h = 20. 
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Type of data 


Method 


Twist angle, deg 




0 

15 

30 

45 

60 

Finite element 

A 

B 

C 

D 

E 

F 

G 

H 

I 

J 

K 

L 

M 

N 

0 

3.399 

3.396 

3.422 

3.287 

3.312 

3.354 

3.417 

3.394 

3.419 

3.385 

3.383 

3.351 

3.349 

3.472 

3.287 

3.310 

3.312 

3.413 

3.318 

3.315 

3.313 

3.311 

4.024 

3.293 

3.320 

3.330 

3.411 

3.403 

3.286 

3.287 

3.375 

3.287 

3.320 

3.416 

3.412 


3.406 




3.394 

3.609 
3.397 
3.384 

3.610 

3.366 

3.287 

3.167 

3.609 

3.394 

3.378 

3.630 

3.020 

3.397 

3.382 

3.620 

3.396 

3.381 

3.620 

3.389 

3.375 

3.640 

Shell 

P 

o 












V 






Beam 

R 

S 

3.342 

3.354 

3.344 

3.357 

3.349 

3.364 

3.357 

3.378 

3.369 

3.396 

Experiment 

T 

U 

3.213 

3.213 

3.192 

3.200 

m 

3.190 

3.198 

3.222 

3.226 

Difference between 
average of the- 
oretical results 
and experimental 
result T, percent 

All 

5.66 

5.90 

5.10 

5.63 

5.01 

Finite element 
methods (A to O) 

5.88 

6.15 

5.20 

5.70 

5.05 

Selected methods 
(D and H to N) 

5.83 

4.94 

4.34 

5.26 

2.20 
































































TABLE III.— Continued. 


(35) First torsion mode: 1/0; a/b = 3; b/h- 5. 


Type of data 

Method 


Twist angle, ^ 

J * deg 




0 

15 

30 

45 

60 

Finite element 

A 

19.800 

19.860 

20.050 

20.340 

20.740 


B 

19.500 

19.560 

19.740 

20.040 

20.430 


Q 

21.220 


21.450 


22.170 






D 

20.700 

20.740 

20.900 

21.150 

21.510 


E 

21.230 

21.400 

21.900 

22.740 

23.750 


F 

21.230 


21.230 


23.250 






G 

H 

21.690 

19.240 

21.510 

21.420 

21.240 

20.990 

20.060 



21.270 


I 

21.360 

21.680 


J 

K 




19.240 

19.310 

19.530 

19.920 

20.410 


L 

18.510 



18.960 








M 

19.470 

19.550 

19.810 

20.200 

20.710 


N 

18.560 

18.610 

18.750 

19.000 

19.340 


O 

19.460 

19.510 

19.680 

19.950 

20.310 

Shell 

p 






Q 












Beam 

R 

19.270 

19.340 

19.560 

19.910 

20.390 


S 

19.090 

19.170 

19.390 

19.780 

20.300 

Experiment 

T 

18.430 

18.450 

18.670 

18.890 

19.320 


U 

18.460 

18.500 

18.700 

18.910 

19.460 

Difference between 
average of the- 

All 

7.73 

7.14 

8.18 

6.80 

8.30 

Finite element 
methods (A to O) 

8.25 

7.78 

8.83 

7.19 

8.79 

oretical results 

and experimental 
result T, percent 

Selected methods 
(D and H to N) 

5.89 

5.64 

6.89 

4.82 

6.30 















































































Type of data 


Method 


Finite element 


Shell 


A 

B 

C 

D 

E 

F 

G 

H 

I 

J 

K 

L 

M 

N 

O 


P 


Beam 


Q 


R 

S 


Experiment 


T 

U 


Difference between 
average of the- 
oretical results 
and experimental 
result T, percent 


All 


Finite element 
methods (A to O) 


Selected methods 
(D and H to N) 


III. — Continued. 
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Type of data 


Method 


Twist angle, deg 






0 

15 

30 

45 

60 

Finite element 

A 

110.200 

110.300 

110.700 

111.200 

112.000 


B 

108.200 

108.400 

108.700 

109.200 

110.000 


C 

121.000 


121.400 


122.600 


D 

115.500 

115.500 

115.600 

115.900 

116.300 


E 







F 







Q 







H 

109.200 





110.500 


I 

j 

120.000 


119.600 


117.600 


K 

105.400 

105.000 

103.500 

100.800 

97.460 


L 

101.500 


— 

102.400 



M 

107.000 

106.900 

106.700 

106.700 

107.700 


N 

99.130 

99.220 

99.480 

99.920 

100.500 


0 

107.000 

108.100 

108.400 

109.000 

109.700 















































TABLE III. — Continued. 


(43) Fourth torsion mode: 1/3; afb~ 3; b/h- 5 
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Type of data 


Method 


Twist angle, <p , deg 


Finite element 


Shell 


Beam 


Experiment 


Difference between 
average of the- 
oretical results 
and experimental 
result T, percent 



A 158.500 158.900 160.100 161.700 163.000 

B 158.500 158.900 159.900 161.300 162.200 


































Frequency parameter 








Twist angle, deg 


(25) Second bending mode: 0/1; a/b = 3; b/h = 20. 

(26) Third bending mode: 0/2; alb = 3; b/h = 20. 

(27) First edgewise bending mode: -/-; a/b = 3; b/h = 20. 
(28) Second torsion mode: 1 / 1 ; a/b = 3\ b/h = 20. 

Figure 4.— Continued. 




Twist angle, deg 


(29) Third torsion mode: 1/2; a/b = 3; b/h -20. 

(30) Fourth bending mode: 0/3; a(b- 3; b/h = 20. 

(31) Fourth torsion mode: 1/3; a/b = 3; b/h = 20. 

(32) Chordwise bending mode: 2/0; a/b = 3; b/h = 2. 

(33) First bending mode: 0/0; a/b = 3; b/h = 5. 

(34) First edgewise bending mode: a/b = l; b/h = 5 

Figure 4.— Continued. 
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Figure 5 — Continued. 
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Method Mode 
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(3) Concluded. 
Figure 5 — Continued. 
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(5) a/b = 3; b/h = 20\ <5 = 0 
Figure 5 — Conlinued. 
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(7) a/6 = 3; b/h = 5; * = ( 
Figure 5 — Continued. 
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